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It is believed that soon after the Planck time, Einstein's general relativity theory should 
be corrected to an effective quadratic theory. In this work we present the 3+1 decompo- 
sition for the zero vorticity case for arbitrary spatially homogenous spaces. We specialize 
for the particular Bianchi / diagonal case. The 3— curvature can be understood as a 
^ ■ generalized potential, and the Bianchi / case is a limiting case where this potential is 

^^ ' negligible to the dynamics. The spirit should be analogous, in some sense to the BKL 

^^, ' solution. In this sense, a better understanding of the Bianchi / case could shed some 

jys , light into the general Bianchi case. 

rr\ • 

, ■ Keywords: Quadratic gravity; effective gravity; homogenous geometries. 
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1. Introduction 



The semi-classical theory consider the back reaction of quantum fields in a classical 
^ , geometric background. It began about forty years ago with De Witt Ref. 1 , and since 

C^ ' then, its consequences and applications are still under research, see for example Ref. 

m 

Different from the usual Einstein-Hilbert action, the one loop effective gravita- 
tional action surmounts to quadratic theories in curvature, see for example Refs. 
[U [H It is the gravitational version of the Heisenberg-Euler electromagnetism. As it 
is well known, vacuum polarization introduces non linear corrections into Maxwell 
electrodynamics, first obtained by Heisenberg- Euler. 

This quadratic gravity was previously studied by Starobinsky.'^ It is of interest 

for example in the context of the final stages of evaporation of black holes, infla- 

[7I • • • ^ 

tionary theories ,'-'-' in the approach to the singularity'-' and also in a more theoretical 

context.'^ 

The effective gravity, was apparently first investigated in Tomita's article Ref. [10] 

for general Bianchi / spaces. They found that the presence of anisotropy contributes 

to the formation of the singularity. Berkin's work shows that a quadratic Weyl 

theory is less stable than a quadratic Riemann scalar R^. In particular, in the very 

interesting article Ref. 11, Barrow and Hervik addressed the anisotropic cases of 
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Bianchi / and //. The most interesting results of Barrow and Hervik are the exact 
solutions for quadratic theories of the same type investigated in this present article. 
Instead of the metric, the field equations in Ref . [TT] are written in a different set of 
variables which by now is a well known procedure used in the context of dynamical 
systems in cosmology!^. Homogenous solutions in the context of quadratic gravity 
was also addressed by Ref. [131 There is an interesting article by Saridakis, Ref. [M] 
in which the Kantowski Sachs anisotropic case is addressed for quadratic gravity. 
Schmidt does a review of higher order gravity theories in connection to cosmology.^'^ 

Also in the context of quadratic theories we have the alternative of Gauss-Bonnet 
type F{G) by Odintsov, Nojiri and collaborators Ref. [TBI (for recent reviews see Ref. 
[Tfl) . When the action depends on a arbitrary function of the Gauss-Bonnet term, it 
is not a top invariant and a consistent dynamic follows from it. Theories of E? type 
are also investigated, for example in Ref. [TH 

In our case the Gauss-Bonnet term is understood as a surface term and dis- 
carded. 

We have previously addressed the Bianchi / solutions in Refs. [El [20l Stability 
is also discussed there. More recently we have submitted two articles, one in which 
we investigate how the approach to Minkowski space occurs in the Bianchi / case. 
l^ In the other article we present the Bianchi VII a solutions Ref. [21] In Ref. [22l 
it is obtained that the solution is a superposition of a pure tensorial component, 
and a pure scalar component. Thus the approach to Minkowski space involves the 
production of scalar and tensorial gravitational waves. Speculations about the va- 
lidity of the semi-classical argument can be raised, anyway, in a certain sense the 
quadratic counterterms which we consi der are the most natural ones expected from 
the renormalization of a quantum field. l^ 

The purpose of this article is to cast the dynamical equations of motion in a 
3-1-1 decomposition for the zero vorticity case. A time like and geodesic vector can 
be defined and we impose the homogeneity on the 3-space. Apparently it is the 
first time the equations are presented in this form. The Bianchi / case is revisited 
but now the equations are written in an analytical fashion. We intend to use the 
expressions obtained in this work in our future works. Most of the article is there 
just for completeness. 

The article is organized as follows. In section 2 a brief exposition of the vacuum 
polarization by the external gravitational field. In section 3 we present the 3-1-1 
decomposition for any Bianchi type. In the section 4 we specialize to the Bianchi I 
case. And the conclusions are presented in section 5. 



2. The Divergent Counterterms 

Since the contribution of other spin fields to the effective action, are of the same 
type for the scalar field, see for example Ref. I23[ we will considered a quantum 
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scalar field in a curved classical background only 

After an integration in parts in the above action, the gaussian integral results as 



e 



= I D4>cxpl^ I d^Xy/^(l)[\/a\/'' + m^ +^R](f>\ 



iW 



[det(iF/27rM2)]-i/2 
W = -ln(det(iF/27rM2)) = -tr \n{F) + C, 

where {x\F\y) = (VoV" + m^ + £_R)S{x — y)/\/—g, M and C are constants, and 
the trace means tri3 = linij,_i.y J d'^Xy''—g{x\B\y). The effective Lagrange function 
is obtained as a Mellin transform of the Schrodinger kernel 

SW^'-tr[^)^ltry^^dse-^-SF 

i ( r°° e"*"-^ 
dW = --Strl / ds 



where the kernel K{s, x) — Imix^xi K{s, x, x'), satisfies a Schrodinger type equation 

d 

OK(s, X, x') + ^RK(s, X, x) + m K(s, x, x) ~ i-—K(s, x, x'), 

OS 

together with the boundary condition linis^o K{s, x, x') — d(x — x'). This equation 

can be solved perturbatively, 

oo 

K{s,x) - lim ir(s,x,x') =---^e-*('^/2«+"'^) Va„(2sr 

^ ' n=0 

where a is half of the length squared of the geodesic 



1 [\ ( dx-dx^y 



and the a„ are known as Seeley-de Witt coefficients. The ones connected to diver- 
gencies are just the first three 

'^abcdu _ _ Dab D _\ ^ I '' c\ Y7aY7 D 1 I C 1 d2 



/ /'OO T 2 



2 
n=0 
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The term Va^'^R is a total derivative and will be omitted. Also making use of the 
Gauss-Bonnet surface term 



/ d^X^^{R'''''''Rabcd - ^R'^'Rab + i?'), 



the square of the Riemann tensor term can be written in a convenient combination 
RabR"'^ ~ l/3i?^ which is dynamically equivalent to a term proportional to the 
square of the Weyl tensor Cabcd.C"'^^'^. A theory without these specific counterterms 
is certainly inconsistent, except for very particular situations. 

With this in mind, the field equations for the semiclassical theory, are obtained 
performing metric variations in the gravitational Lagrangian 



C = ^/^ 



-k + R + a[ RabR"^ ~ ^^' 1 + ^^' 



(1) 



where a and /? are constants. Classical sources are not taken into account in this 
present work. For the spatially homogenous space they are described by the tensor 

E = EabUj" ® UJ^, uj° = dt, 

Eab - Gab + IgabA - ((3 ~ ^o) F^ - aH^^^ = 0, (2) 

where 

Gab — Rab — -^QabR, 

^ab = '^9abR — 2RRab — ^Qab^R + "^R^ab, 

^ab — '^9abR'' Red — ^Rab — -^Qab^R + R-ab — 2i?'^ Rcbda- 

First let us emphasize that every Einstein space satisfying Rab = Qab^/'^ is an exact 
solution of ([2]), apparently first realized in Ref. [TTl As any other metric theory, 
the covariant divergence of the above tensors are zero. And the the zero divergence 
S/aE^ = implies that Eqq and £'oq — Oj are in fact constraints, while the dynamical 
equations are contained in the spatial part Ea/s = 0, see for example Ref. i24j p. 165. 



3. 3+1 Decomposition 

We briefly develop the general homogenous case, for more details see Ref. [25l The 
right-invariant basis are the Killing vectors, and satisfy the Lie algebra 

where the greek indices run from 1 to 3. On the other hand the left-invariant basis 
vectors satisfy 

[ea,e/3] = C^/je^ 
[eo,ea] = 0. 
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The connection is defined as 

VaGb = rg„ec. 

The anti symmetric part of the connection follows from 

VfSCa - VaSfj = [e/3, £„], 

which can be regarded as the definition of zero torsion. Metricity Vcgah ~ implies 
that 

Tatic = -Z \9ab\c + 9ac\h ~ 9hc\aj + „ \~^ahc + (-^hac + ^ cah) i 

where the | means the directional derivative /|a = e\di] . We restrict the analysis 
to rotation free models in this present work 

where w" is the left invariant 1— form basis, dual to e^. Once the structure constants 
are known, in fact, the 1— form basis can be found by solving Cartan first structure 
equation 

The time like vector field u" — (1,0,0,0) is geodesic and rotation free, as can be 
easily checked for the above line element and connection. In this particular case the 
extrinsic curvature 

K-afi — '^aU[j — T^fj — -ha[3 
i\^ — II J:\^fia — ^ aO — ^ Oa^ 

and the spatial part of the connection reads 
The 3 and 4 Riemann curvature follows 

3 pa T^a T^P j^a -pP /^p j^a 

-"■/3/ii/ — ^ piJ.^ I3u ^ pv^ /3^ '^pv^ Pp 

R'pnv = ^R'ppv + K'^Kfiy - K"Kj3f, 

^al3p = ~^f3p^ap + KppV^^ - K p^T ^^^ 

note that K — dt {K\ + Xl + Kf). The Ricci tensor and Riemann scalar follow 
from contractions of the above tensors 



a/3 


= ""R^p + KK^p - 2KP„Kpp 


Oa 




00 


= -k-K,jK'^ 


= 


^R + 2k + K,jK^^ + K'^. 



K, 



a/3 
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This theory depends on the derivatives of the Ricci tensor and Riemann scalar 

DR= -R-KR 

DRoo = -i?oo - RooK + 2R^pK'^KPf + 2K'^C$^R^ + 2C^^7^"''i?o/3 

+2K,^X"^i?oo 

Di^oa = -iJoa + '^R^pK - KRqo. + 2if '"^r^^i?,^ + i?00 (-C^.^p" + ^^p^a) 

+R„^K^''K^, + i?o. [K + KK^ + K^r-^y^ + C;^KphP^) + 2Ro.K'''^K^^ 

aR^i3 = -^Rc.p + 2RpfiK - ^KR^fi + i?.prs^r^^/i^T + RooKk^^ 
-R.^Ki^K^p + Rp, (x^ + KK'^^ + rs^r;;^/^^^ + c>;^Koh'' 



after the substitution of the above expressions into ([2]), it can be seen that the 
higher time derivatives are cointained in the spatial parts of H^^^ and H^^\ 



4. Bianchi I 



reason is that the ^Rp^g, acts as a generalized potential much in the same sense 
as in the mixmaster case. In this sense, ^Rg g, can decrease and become irrelevant 



We shall apply the preceding expressions to the particular spatially flat case. The 

and the dynamics can be arbitrarily close to the Bianchi / case. 

This is the abelian case in which all the structure constants are zero. Of 
course '^Rap^.v = 0. We will restrict ourselves to diagonal extrinsic curvature 
Kap = diag [Ki,K2tK^], which is consistent with a diagonal metric also hap = 
diag [af , aj, a\]. Note that Kap.h'^'^ = S^{Hp + 2(_ff^)^), where H/j = ap/ap are the 
Hubble constants in each direction, and as usual, the trace K — Hi+ H2 + H3 ~ 9, 
is the expansion 



Ri=5i(eHp+Hp^ 

Roa = 

R = 29 + e^+H„Hp6°''^. 
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In the following we quote the result for the contribution coming from the R^ term, 



-^00 - ^2 '^^ 



H^HJ^'' +2{2e + e' + H^HJ'^'' {6 + H^.H^S^^") 



H 



40{9 + 60 + Hf.HJ"") 



H^^^ -0 



(1)" 



^2 



29 



HuH^Sf'n -2 26l + r + i/ui?, 



W"") (■ 



9Ha + Ha 



+4 



{0 - Ha) (6 



H^HJ^'" ]+4 'e 



H^HJt''' + H^H.S^'^ 



the spatial trace of which, i/*^L" = H^^^^ + H^'^^^ + H^^^s , 



i/Ci)^" ^^2^ + 612 + H^HJ^"A + 129 + 89^ + 2099 + 29^9 - 29^ 
-2Hf,H^S'"' (9^ +9)+ WH^HJf" + UH^H^dt"" + 12H^HJfr 

And the contribution from the RabR"^^ term 



H 



2 - Hf.HJ^^'' 



d + H^H^S^^'j +J2[0Hr. + H 
H^.HJ^" + 2H^,HJ^'' {9 + HpHJP^^ + 2H^ (9H, + ij,) 5^''' 



SHpHJP''^ 



H^2 - 



H 



(2)« 



5% 



1 



9 + H^H,,S'"') +^{9H^ + H^ 



2{Haf {9 + H^H.S^"' ] + {9Ha + 29Ha + 9Ha + H, 



2Ha [9 + 99 + Hf.HJ'"' 



-2 



+9 [9 + 99 + Hf.HJ^"' 
+ {9 + HpH.S'"'^ {ka + (H 
with spatial trace 



/7(2)a ^ £ 

"2 



h 9^ + 99 + Hf.HJ'"' + H,,H,S^"' 
HJ^^'^H, {9H^ + ij^) 

h 9 {9Ha + 9Ha +H^Y 



769 + 39^9 



9 + H^HJ'^'^f + J2 \0H, + H,^ 
+3^2 _^ 4-g,- _^ ■iH^HJf''' - 9H^HJ>''' - 2H^HJ^"' (9 + 9^ + SHf.H^ud'"' . 
Now we can use the sum of the spatial part of field equations (HJ to obtain the higher 
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derivative of the expansion 9, giving the equivalent of Raychaudhuri's equation 

"e = — {-29 - -9^ - -H^HJ"" + -A 
12/3 \ 2 2^2 



((3-^) ^ h9 + 9^ + H^.HJ^'A + 8^2 + 20616' + 29'^9 - 29^ 
2Hf,HJ'"' (9^ + 9)+ 89Hf,HJ'"' + UH^HJ^"" + UHf.HJ'"' 

J2 {SH^ + H^) ) + 799 + W'^9 



H^HJ^'^f 



-36*2 + iHaHJ^"" - 9Hf,HJ'"' - 2i/„iJ, 



\,s'"' (9 + 9^) 



ZH.H,,u5^''' 



Again, using the field equations ^ the higher derivatives Ha can be isolated, and 
9 can be substituted in the following expression yielding a consistent dynamical 
system 

H„ = i l9Ha + Ha-9-^- iff^H.5^'' + ^A 



^ 3 



1(2.- 



-^{6 -Ha) 



9' + H^H^Sf^'j -2 [29 + 9' + Hf^HJ^'^j [9Ha + H^ 
Hf.HJ^'') + 4 ('9 +9^ + 99 + H^.HJ^" + 2H^HJ'"' 

^2{Haf(9 + H^,HJ^''' 



d + H^H^S^^n +^(9H, + H, 



1 
2 

+ ( 9Ha + 29Ha + 9Ha ] - 2Ha (9 + 99 + H^HJ'"' 



9^ + 99 + Hf.HJ^'" + Hf.HJt"' 



'H, 



J^^H, ( 



9Hi^i + H^ 



-2 

9Hn + 9Hn + H, 



H^HJ^^'^ 



H^HJ'^'^ ]{Ha + [Ha 



)]} 



In the above expressions only the spatial part of ([2]) is used, while not shown here is 
the 00 component of Eat which acts as a hamiltonian constraint. Let us emphasize 
that 9 — Hi + H2 + H3, thus this last two equations are not independent. 



5. Conclusions 

The semi-classical theory consider the back reaction of quantum fields in a classical 
geometric background. It began about forty years ago with De Witt Ref.[Tl and since 
then, its consequences and applications are still under research, see for example Ref. 

m 

Different from the usual Einstein-Hilbert action, the one loop effective gravita- 
tional action surmounts to quadratic theories in curvature, see for example Refs. 
[U [H It is the gravitational version of the Heisenberg-Euler electromagnetism. As it 
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is well known, vacuum polarization introduces non linear corrections into Maxwell 
electrodynamics, first obtained by Heisenberg- Euler. 

In this article we present the decomposition 3 + 1 for arbitrary spatially homoge- 
nous space times for this particular effective quadratic theory. It turns out that for 
this particular gravitational theory ([1]), every Einstein space satisfying Rab = gab-^/'^ 
is an exact solution, this includes of course the vacuum case A = 0. Anyway the 
most interesting is the BKL oscillatory approach to the singularity which occurs for 
vacuum in Einstein's context !^^ Since the mixmaster solution is a vacuum Rab = 
solution for the Bianchi IX case it will occur exactly the same in the effective theory 
considered in ([T|). 

In this sense a better understanding of the Bianchi / case for the quadratic 
theories seems attractive. The 3 curvature also can be understood as a potential 
and while approaching the singularity its influence in the dynamics can decrease 
almost to zero. Thus, an arbitrary Bianchi solution of the quadratic gravity can 
approach very much the Bianchi / case. We intend to use this 3 + 1 decomposition, 
numerically, in future works. Particular Bianchi / analytic exact solutions for a 
quadratic theory identical to this one were already found by Ref. [TTJ Also for f{R) 
gravity models including the R + R^ case, Bianchi I space-times were studied in 
Ref. 1271 where it was shown that equations for the anisotropic part of the metric 
can be integrated. 

We emphasize that the approach of this present work is facing the theory ([T]) as 
an effective, and classical theory. While considering it as a candidate for a quantized 
gravity does introduces ghosts. Which together with tachyons and the additional 
degrees of freedom in contrast to Einstein's theory, makes it improbable that the 
BKL solution will be the generic one in this case. 
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